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Abstract: We study the localization length of few interacting particles in a random 
potential. Concentrating on the case of three particles we show that their localization 
length is strongly enhanced comparing to the enhancement for two interacting particles. 

PACS. 72.15Rn - Quantum localization 
PACS. 71.30+h - Metal- insulator transition 
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Recently it had been shown that in a random potential two repulsing/attracting par- 
ticles can propagate coherently on a distance l c which is much larger than one particle 
localization length l\ in absence of interaction |ffl]. In some sense interaction destroys 
quantum interference which leads to one-particle localization and creates an effective pair 
of two particles of size l\ propagating on a large distance. For better understanding 
of this result Imry developped [Q] a scaling block picture of localization for interacting 
particles which can be applied in principle for a larger number of particles and higher di- 
mensions. Intensive numerical investigations of Pichard and coworkers [|| and von Oppen 
and coworkers Q confirmed the existence of the two interacting particles (TIP) effect. 
While some additional checks are still required the results H [|| definitely show that in 
one-dimensional case the TIP length is l c oc l\ a with a close to the theoretically pre- 
dicted power (lj a = 2. These results are also in agreement with the previous studies of 
Dorokhov who analysed the case of two particles confined by strong attraction in a well 
with size much smaller than l\ ||. The investigations of TIP effect in higher dimensions 
have been done in B, ||, [5J and they demonstrated that in dimension d = 3 the TIP pair 
can be delocalized below one-particle Anderson transition where all one-particle states 
are localized. 

While now the properties of TIP propagation reached a level of qualitative understand- 
ing, the problem of a larger number of interacting particles is still not well understood. 
From the physical point of view the most interesting situation is the case of finite density 
of particles. However, the analysis in this case is quite complicated and at present only 
estimate [^| and numerical studies in Q have addressed this problem. One of the ways to 
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approach this problem is to analyse the case of larger number of particles. The simplest 
case is three interacting particles where the situation is not so trivial since two particle 
interaction leads to the Breit-Wigner structure of eigenstates || [To|, 11, 12]. In this paper 
we will concentrate on this three particle model. 

However, before the analysis of three particle model let us first discuss some general- 
ized model of TIP (see also []13| ). In this model the first particle is moving on sites n\ 
in a one channel Anderson model with diagonal disorder changing in the interval ±Wi, 
the intersites hopping matrix element is V\ and localization length at the center of the 
band is 1% ~ 25(Vi/Wi) 2 1. The second particle is moving in a strip with M transverse 
channels with sites marked by index ni along the strip and index ^(l < ^2 < M) in 
transverse direction. The disorder in the strip is independent of disorder in the chain with 
first particle and the localization length for the second particle is 1% oc M. The hopping in 
the strip is V ~ V±, and we assume that 1% > l\. Now we will analyse what will happen if 
the interaction of the form U5 ni ^ n2 is switched on between two particles. Similarly to Q 
one should first estimate the transition matrix elements U s between eigenstates without 
interaction ( U = 0). This gives 

Us = U ^ ^ni,mi^n 2 ,n2,m2,m2^ni,m' 1 ^n2,n2,m2,m' 2 ^ n i' n 2 W 
ni,n 2 ,n 2 

where R represents the transformation between the lattice basis and one-particle eigen- 
states so that Rn urni exp(-| n\ - mi \/h - ^m.mj/v^l and Rn 2l h 2 ,m2,rh2 ~ 
exp(— | 722 — Tt2 \/h — i6n 2 ,n 2 ,m 2 ,m 2 ) / v ' MI2 correspondingly for the first and second par- 
ticle. The phase 9 randomly changes with indices. Due to the exponential decrease of R 
one should take into account only the states with | n\ 2 — m\ 2 |< l\2- For the case I2 > h 



the sum in (|]) containes approximately l\M random terms so that U s ~ U / foy/hM). 
The interaction induced transition rate is given by the Fermi golden rule T ~ U s 2 p c where 
Pc ~ l\l%MjV is the density of coupled states. As the result, T ~ U 2 /(VI2) is independent 
of h and M. 

With the rate V we can determine the interaction induced diffusion rate for the first 
particle which is D\ ~ l\ 2 T ~ U 2 l 2 /{Vl 2 ) and appears as the result of collisions of the 
first particle with the second one oscillating in the block of size ^2- Knowing the diffusion 
rate, it is possible to determine the localization length for a pair in a way similar to 
that used for the kicked rotator [|14|] and based on the uncertainty relation between the 
frequency and time (see also Indeed, the number of excited states in the first chain 

grows with time t as An\ ~ (D\t) l l 2 . Since two particles are propagating together so that 
I n\ — n% |< I2 the total number of excited states in both chains is AiV ~ An\(Ml2)8E /V 
where 5E takes into account the factor that the states are excited only in some energy 
interval inside the band width V. Generally, 5E < V and it is of the order of Breit- 
Wigner width T |J, but we will see that 8E does not enter in the final expression for 
the localization length of pair, and therefore actual value of 5E is not very important 
(see also Indeed, all these A./V levels are homogeneously distributed in the energy 

interval 5E and the average splitting between them is Av ~ 5E/AN. According to 
the uncertainty relation between frequency and time at the moment t we can resolve 
discrete lines with the splitting 1/t. Therefore, at the moment t* defined by the equation 
Au ~ 1/t* the discreteness of the spectrum is resolved and the diffusive propagation is 
stopped at t* ~ AN(t*)/5E. This condition gives the localization time t* for TIP pair 
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and the localization length for the first particle l c \\ 

t* ~ U 2 h 2 M 2 l 2 /V; l c i ~ Am ~ {U/V) 2 h 2 M (2) 

The interesting feature of this result is that l c \ is independent of l 2 . This indicates 
that the nature of motion in the second chain does not influence much the localization in 
the first chain. For M = 1 the length l c \ is the same as in the case of TIP localization 
in one chain. However, the growth of the number of channels M in the strip leads to the 
increase of l c \. The localization length for the second particle is l c2 ~ h if h 3> lci and 
lc2 ~ Id if h *C l c i- The similar approach can be used for analysis of TIP localization in 
higher dimensions |Oj. 

Let us now consider three interacting particles in one-dimensional chain with on site 
interaction Z7i2^m,ri2) ^23^n 2 ,n3 an d Ui35 nijTl3 where m,2,3 marks the site position of cor- 
responding particle in the chain. As above, the one-particle localization length is l\ and 
the band width is 4V. For simplicity we will assume that U13 = and U23 > U\i- Then 
in first approximation the particles 2-3 form a pair of size l\ which is localized on the 
length l c2 ~ {U 2 z/V) 2 h 2 . When this pair approaches the first particle at a distance l\ 
the interaction between three particles in a block of size l\ gives mixing between /i 3 3- 
particle states. An effective matrix element U s \ of interaction between 3-particle states in 
the block of size l\ should been calculated in the second order perturbation theory, since 
direct interaction couples only 2-particle states. Therefore, the matrix element between 
initial state 1 123 > and final state |1'2'3' > is given by diagram presented in Fig.l with 

intermediate state |1'23 >. It is of the form 

u _ v <12\U 12 \1'2><23\U 23 \2'3' > _ U 12 U 23 
sl ^(E 1 +E 2 + E 3 -E v -E 2 -E 3 ) ~ l 1 3 A 1 U 



It is important that the summation is carried out only over single particle states 2, hence 
Ax ~ V/h is single particle level spacing. Finally this gives the mixing rate in a block of 
size l\ 

T l ~ U sl 2 p 3 ~ (U 12 U 23 /V 2 ) 2 V/h (4) 

where p 3 ~ /V is the density of 3-particle states in the block. This Ti gives the 
mixing rate during the collision of the first particle with the pair 2-3 in the block l\. 
The frequency of such collisions is of the order of h/l C 2 since from ergodicity the ratio of 
time of the collision to the time between collisions is proportional to the ratio of volumes. 
Therefore, the average transition rate for 1-particle per unit time is Ti ~ Tili/l c2 . Such 
transitions give the diffusion rate of the first particle D\ ~ T\li 2 ~ U\ 2 2 /V since the size 
of transition is l\. Similarly to the previous case with two chains the total number of 
excited states after time t* is AiV ~ (Dit*) l / 2 {l c2 h)5E /V where SE is an energy width 
in which the levels are mixed. The localization time t*, as previously, is determined from 
the condition AN ~ 5Et* which gives 

^Oi(Wi) 2 /^ 2 ; l*/h ~ D x l c2 /V ~ (U 12 U 2 3/V 2 ) 2 l 1 2 (5) 

For Ui 2 ~ U 2 3 ~ U the localization length for the first particle is enhanced only if 
there is an enhancement for two-particle localization length, namely {U /V) 2 l\ > 1. This 
result is quite natural since for (U /V) 2 l\ < 1 two-particle interaction is too weak and 
it is not able to mix three-particle levels. Another limiting case in (||) corresponds to 
U ~ V. For such interaction l c \ ~ Zi 3 which is similar to the case of three particles 
trapped in a bag of size l\. Indeed, one can consider 3-particle bag model like TIP one 
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HI with effective number of transverse channels M e f ~ l\, therefore for the 3-particle bag 
/fe3 ~ M e fl\ 2 ~ The same estimate for ^3 had been also obtained in basing on the 
approachs developped in || |2|. In some sense the result @ shows that similar to the TIP 
case the "size" and "form" of the bag is not important for the effect. Let us also mention 
that the case U23 ~ V is similar to previously analysed model of TIP in the chain and the 
strip (^). Indeed, here the third particle gives the effective number of channels M ~ l\ so 
that (|5|) becomes equivalent to (Q). Generalization of the result (||) for k particles gives 
the enhancement l c k/h ~ ({U/V) 2 h) k ~ 1 . 

For the 3-dimensional case l\ in the enhancement factor (U /V) 2 l\ should be replaced 
by Z1 3 @, |, m so the derealization takes place if ((U/V) 2 h 3 ) k - 1 > 1. This means that 
the derealization border for few particles coincides approximately with that for TIP and 
therefore it is not possible to have propagating cluster with k > 2 repulsive particles. In 
some sense only TIP pairs are well defined. 

We thank J.-L.Pichard and Y.Imry for useful and stimulating discussions and com- 
munication about their results prior to publication. One of us (DLS) gratefully 
acknowledges the hospitality and fine working conditions at the Institut Henri Poincare 
during research on the above problem. Another coauther (OPS) thanks Laboratoire de 
Physique Quantique, Universite Paul Sabatier for hospitality and financial suport during 
the work on this problem. 
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Figure captions 



. 1: Diagram for effective three particle matrix element U s \ 
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Fig. 1 



